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' Abstract 

| We prove edge universality of local eigenvalue statistics for orthogonal invariant ma- 

trix models with real analytic potentials and one interval limiting spectrum. Our start- 
ing point is the result of [21j on the representation of the reproducing matrix kernels 
of orthogonal ensembles in terms of scalar reproducing kernel of corresponding unitary 
f-^ | ensemble. 

1 Introduction and main results 

We study ensembles of n x n real symmetric (or Hermitian) matrices M with the probability 
distribution 

P n (M)dM = Z-J, exp{-^Try (M)biM, (1.1) 

where Z nt p is a normalization constant, V : M — > M+ is a Holder function satisfying the 
^ . condition 

(N 



\V(X)\ >2(l + e)log(l + |A|). (1.2) 

A positive parameter (3 here assumes the values /3 = 1 (in the case of real symmetric matrices) 
or (5 = 2 (in the Hermitian case), and dM means the Lebesgue measure on the algebraically 
independent entries of M. 

The joint eigenvalue distribution corresponding to (|1.1|) has the form (see 



OO 



PnA^-An) = Q^U e ~ nmXi)/2 II (1-3) 
i=l l<j<fe<« 

where Q ny p is a normalization constant. For both cases (j3 = 1, 2) the behavior of Normalized 
Counting Measure (NCM) of eigenvalues is now well understood. According to NCM 
converges weakly in probability to the non random limiting measure known as Integrated 
Density of States (IDS) of the ensemble. The IDS is absolutely continuous, if V satisfies the 
Lipshitz condition [19]. The non-negative density p(A) is called Density of States (DOS) of 
the ensemble. IDS can be found as a unique solution of a certain variational problem (see 

EHISl). 

To study the local regimes for ensembles (jl.lj) means to study the behavior of marginal 
densities 

p\ n )(Xi,...,Xi) = p ni/3 (Xi,...Xi,Xi +1 ,...,X n )dXi +1 ...dX n (1.4) 

in the scaling limit, when Aj = Xo + Xi/n K (i = 1, . . . , I), and k is a constant, depending on the 
behavior of DOS p{X) in a small neighborhood of Aq- If p(Ao) ^ 0, then k = 1, if p(Xq) = 
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and p(X) ~ |A — Xo\ a , then k = 1/(1 + a). The universality conjecture states that the scaling 
limits of all marginal densities are universal, i.e. do not depend on V and depend only on a 
and j3. One of the most known quantity probing the local regime is the gap probability, i.e., 
the probability that there is no eigenvalues in the interval A n (a, b) = [Ao + a/n K , Ao + b/n K ] 

E n ^(A n (a,b)) = e{ H(l - l An(a ,6)(A t ))}. (1.5) 
k=l ) 

Thus, results on the universality of local eigenvalue statistics usually include proofs of uni- 
versality of the gap probability. 

For unitary ensembles all marginal densities can be represented (see [H])) in terms of so 
called reproducing kernel 

iUA,M) = EW n) m (n V)- (1.6) 

1=0 

where 

^ (n) (A) = e W {-nV(X)/2}p ( l n) (X), I = 0,..., (1.7) 

and {p[ n ^}f =0 are orthogonal polynomials on R associated with the weight w n (X) = e~ nV (^) ) 
i.e., 

P < i n) (X)p^(X)w n (X)dX = 5 l , m . (U 



In particular, 

E n>2 (A n (a,b)) = det{l - #A„(a,6)}, 

where det{. . . } is the Fredholm determinant and K A r a ^ is the integral operator with the 
kernel (jl.6l) in L 2 (A n (a, b)). Hence, the problem to study marginal distributions is replaced 
by the problem to study the reproducing kernel K n (X, fx) in the scaling limit. 

The problem was solved in many cases. For example, in the bulk case (p(Xq) ^ 0) it was 
shown in [15] (see also [T7]) that for a general class of V (the third derivative is bounded 
in the some neighborhood of Ao) the scaled reproducing kernel converges uniformly to the 
sin-kernel. This result for the case of real analytic V was obtained also in [7]. 

Universality near the edge, i.e., the case when Ao is the edge point of the spectrum and 
p(X) ~ |A — Aoj 1 / 2 , as A ~ Ao, was studied in [7]. It was proved that 

lim ) K n 2 (X + si/n 2/3 7, A + s 2 /n 2/3 j) = Qai(si, s 2 ), 

where 

Ai'(si)Ai(s 2 ) - Ai( Sl )Ai'(s 2 ) 

QAi{si,s 2 ) = . (1.9) 

si - s 2 

This result for GUE (V(A) = A 2 /2) was obtained in [26J. There are also results on universality 
near the extreme point, where p(X) ~ (A — Ao) 2 , as A ~ Ao (see [1] for real analytic V and 
[20] for general V). 

For orthogonal ensembles {(3 = 1) the situation is more complicated. Instead of (jl.6p we 
need to use the matrix kernel 

K ^ X ^ ) -\lS n {X,p)-e{X-p) S n (p,X) )■ (L10) 

Here 

n-l 

S n (X,p) = -J2 4 n \*)(M (0 ' n) );](ne^ n) )(p), (1.11) 

i,j=0 
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where ip!- are defined by (|1.7p - (|1.8|) and the matrix .A/^ '™) is defined as 

M u = n(^ n \e4 n) ); = {M^ =0 ; = {M^}™^, (1.12) 

where e is the integral operator with the kernel 

e(A) = ^sign(A); e/(A) = J <X- (1-13) 

The symbol d in (jl.lOp denotes the differentiating with respect to fi, and IS n (X,fJ,) means 
the composition of operators e and S n . Similarly to the hermitian case all marginal densities 
can be expressed in terms of the kernel K n (see |25j). In particular, the gap probability has 
the form 

E nA (A n (a, b)) = detail - K n (A n (a, &))), (1.14) 

where K n (A n (a, b)) is an integral operator from L 2 (A n (a, b)) © L 2 (A n (a, b)) to itself defined 
by the matrix kernel (jl.lOp and det means it Fredholm determinant. The matrix kernel (|1.10|) 
was introduced first in |llj for circular ensemble and then in |14j for orthogonal ensembles. 
The scalar kernels of (jl.lOp could be defined in principle in terms of any family of polynomials 
complete in L 2 (^, w n ) (see |25j), but usually the families of skew orthogonal polynomials were 
used (see [14J and references therein). Unfortunately, using of skew orthogonal polynomials 
for general V rises serious technical difficulties. 

The main technical obstacle to study the kernel (II. lip defined in terms of orthogonal 
polynomials is that there is no uniform bound for ||( J M^°' n ^) _1 ||. According to Widom (see 
[27] ). if the potential V is a rational function, then to control (A4^ 0,n ^) _1 it is enough to 
control the inverse of some matrix of fixed size depending of V (e.g. if V is polynomial of 
degree 2m, then we should control some (2m — 1) x (2m — 1) matrix). In the paper [8] by 
constructing of the exact expressions for the entries of the Widom matrix, it was shown that 
it is invertible in the case V(A) = A 2m . This allowed to prove bulk universality for the case 
V(X) = A 2m + n~ x l 2m ai m -\ A 2m_1 + . . . (in our notations). The same approach was used 
in [9] to prove edge universality and in [10] to prove bulk and edges universality (including 
the case of hard edge) for the Laguerre type ensembles with monomial V. In the papers 
[22\ [23] universality in the bulk and near the edges were studied for V being an even quatric 
polynomial. The most general result for the moment was obtained in [21], where it was 
shown that in the one interval case the matrix, which we need to control, is of rank one. 
This allowed to study any real analytical potential with one interval support and to simplify 
considerably the proof. In the present paper we will use the result of [21J to prove that up 
to some small terms (which do not contribute in the limit) the kernel S is the same that for 
GOE case (see Lemma [1]). This allows us to use the method of [26] to prove universality of 
limiting kernels near the edges. The only difference with [26J is that we use asymptotic of 
orthogonal polynomials of [7], instead of classical asymptotic of Hermite polynomials. 

Let us state our main conditions. 

CI. The support a of IDS of the ensemble consists of a single interval: 

a = [-2,2]. 

C2. V{z) satisfies and is an even analytic function in 

n[di, d 2 ] = {z : -2 - d x < Htz < 2 + di, \53z\ < d 2 }, d 1 ,d 2 > 0. (1.15) 
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C3. DOS p(X) is strictly positive in the internal points A G (—2,2) and p(X) ~ |A =F 2I 1 ' 2 , 
as A ~ ±2. 

C4. T/ie function 

u(X) = 2 I log \n - X\p{p)dp, - V(X) 



achieves its maximum if and only if X G a. 

Note (see |2j) that under conditions C1-C4 the limiting density of states (DOS) p has 
the form 

P(A) = ^P(A)V4-A2 1 |A | <2 , (1.16) 
where the function P can be represented in the form 

= }_ r V'(z)-VVcosy) 

If V is a polynomial of 2mth degree, then it is evident that P{z) is a polynomial of (2m — 2)th 
degree, and conditions C3 guarantee that 

\P(z)\<C, zen[d 1 /2,d 2 /2], P(X)>5>0, AG [-2,2]. (1.18) 

We will use also that under conditions C1-C4 the entries of semi infinite Jacoby matrix 
J^ n \ generated by the recursion relations for orthogonal polynomials (jl.8p 



^<{(A)+grvr^A)+jrvri(A)=A^r(A), jr=o, z= ,.... (i.i 9 ) 

satisfy the relations (see [D [2] ) : qf 1 ^ = and 



r(n) _ -, 

J n+k L 



k 



2nP(2) 



< C ^ 2+ f\ \k\<2n 1 /\ (1.20) 



where P is defined by (|1.17p . Here and everywhere below we denote by C, Co, Ci, c, ... positive 
n-independent constants (different in different formulas). 
The main result of the paper is 

Theorem 1 Consider an orthogonal ensemble of random matrices il.3\) with (3 = 1, and V 
satisfying conditions C1-C4- Set 7 = P 2 / 3 (2), where P is defined by fli. 17\ ). Then for even 
n we have: 

(i) ifp[i is the Ith marginal of il.3\) . then n^^p^ {2+xi/^n 2 ^ , . . . , 2+x//7n 2 / 3 ) converges 
uniformly in Xj > s > —00, j = 1, to the limits coinciding with that for GOE and 
given in terms of 

Q Ai (x, y) = lim n- 2 / 3 7" 1 K n (2 + x/ 7 n 2 / 3 , 2 + y/ 7 n 2/3 ), (1.21) 

n^oo 
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with 

SAi{x,y) 



Q A i(x,y) + \Ai(x) (l-J Ai{z)dz^j , (1.23) 



D A i(x,y) = -d y QAi(x,y) - -Ai(x)Ai(y), 



1 

2' 

Ivi-r-.u) - I QAi(z,y)dy 

Ai(z)dz + / Ai(z)dz I Ai(z)dz 



and QAi(x,y) of fO| 



(ii) if E n x is the gap probability (7777]) °f (E3p> corresponding to the semi-infinite interval 
(2 + s/7n 2 / 3 , oo), then 

Jim^i ((2 + s/ 7 n 2 / 3 ,oo)) := E[ ed9e \s) = det 1 / 2 ^ - Q Ai (s)), (1.24) 

where Q Ai{s) is the integral operator, defined in L 2 (s, 00; w) ©L 2 (s, 00; w^ 1 ) by the 2x2 
matrix kernel HI. 22]) . il.23\) with w(x) = x 2 + 1. 

Remark 1 According to the results of f^j and |i6] /. if we restrict the integration in M.'J\) by 
\W < L = 2 A- di/2, consider the polynomials {pj^'^lfcLn orthogonal on the interval [—L,L] 
with the weight e~ nV and set ipu 1 '^ = e~ nV l 2 p^ ,] "\ then for k < n(l + e) with some e > 

sup |A| < L |^ n ' L) (A) - ^ n) (A)| < e~ nC , \^\±L)\ < e~ nC (1.25) 

with some absolute C. Therefore from the very beginning we can take all integrals in \1.3\) . 
tl.8\) . \1.13\) and h!.12\) over the interval [—L,L]. Note also that since V is an analytic 
function in J) [cf 1,0/2] (see U.15\) ). for any m S N there exists a polynomial V m of the (2m)th 
degree such that 

\V m (z)\ < C , \V(z)-V m (z)\ < e~ Cm , z e Q[d 1 /2,d 2 /2]. (1.26) 

Take 

m=[log 2 n] (1.27) 

and consider the system of polynomials {p(™' L,m ) }^L Q orthogonal in the interval [—L,L] with 
respect to the weight e ~ nV ™^ . Set ^ n,L,m) = p ^' L ' m) e ~nV m /2 and construct ^(M by 
with ip^ 1 ' 1 "' 171 ^. Then for any k <n + 2n^l 2 and uniformly in X £ [-L, L] 



|V>i n ' L) (A) - ^ n ' L ' m) (A)| < e -^i°s 2 -, |#t n ' L) (A) - ^i n ' i,m) (A)| < e~ cx ^ n 
\\M$ n) - M^\\ < e- cl °z 2n , WiM^)- 1 - (X^))" 1 !! < e - c ^ n . 



1.28) 



The proof of the first bound here is identical to the proof of U.25\) (see ]1Q$). The second 
bound follows from the first one because the operator e : L2[—L, L] — ► C[—L, L] is bounded by 
L. The third bound in lil.28\) follows from the first, and the last bound follows from the third 
one and from the fact that \\(Mm' n ^ is uniformly bounded (see f21\j). Hence 

|5 n , m (A, M ) - #n(A,M)l < Cn 4 e- C1 °g 2 - < e - c ' lo ^ n , (1.29) 

and below we will study S n ^ m (X,fi) instead of S n (\,fi). To simplify notations we omit the in- 
dexes m, L, but keep the dependence on m in the estimates. For more detail of the replacement 
see rUF. 
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Our starting point is the representation of the matrix (.M^ '™)) -1 valid under conditions 
C1-C4 (see Corollary 1 in |21j). To formulate this result we introduce a few Toeplitz matrices. 
Consider the infinite matrix V = {.Pj,kYj 3 k=—ao ^ n 00, 00] with entries 

P h k = ^f P ( 2 cos y)e l(j - k)y dy, (1.30) 

and U = V- 1 , 

(On) -1 1 r e i( j~ k } x dx 

'rS > n > = {Rj,k} n jk= Q Rj,k = Rj-k = 7, / stt, c * (1-31) 

J ■>'*■ u J J 2-ir J_ w P(2 cos x) 

It is important that (see [21J), Proposition 1) 

\Rj,k\ < e- c|j - fc| , \((n^ n %l\ < e-^-*!, (1.32) 

where c > is some n-independent constant. Remark also that if we denote by J* an infinite 
Jacobi matrix with constant coefficients 

J* = {Jj t k}jjc=-ooi Jj,k = $j+i,k + fy-i,k, (1.33) 

then the spectral theorem yields that V = P(J*), 71 = P~ l (J*). 
Two more matrices which we use below have the form 

V {0 ' n) = {Dj, k }^l , D hk = 6 j+1 , k - Sj- ltk . (1.34) 
and V(°'°°) = {Vjj}°° l=0 , where 

Here 0(e _c,log2n ) appears because of the integration by parts and bounds (jl.25p . f[1.28]> . 
According to Corollary 1 from [21], under conditions C1-C4 

(M^)~l = Qff + \a,b k + 0(n-V2 i og 6 n) , (L36) 

where 

Q (o,n) = 1 f Vjj 00 ^ for < j < n - 2m, < k < n, Q 3? 

j ' k 2 I (7& ' n >)- 1 X>(°' n )) i)Jfe , for n-2m<j<n, < A; < n. 

0i = ((^°^)- 1 e n _i)„ 6 fc = ((^°>"))-V) fe , r;_ i= 12i (1.38) 

with i?j defined by (OP . 

Note that since (R-)j k = = for | j — k\ > 2m — 2, the standard linear algebra yields 
that (7£(° ,n )) -1 possesses the same property, i.e., 

(n^)~ l k = 0, for \j - k\ > 2m - 2 Q,f^ ] = 0, for |j - k\ > 2m - 2. (1.39) 



and 
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2 Proof of Theorem QJ 

Remark 2 fij is easy to see that the integral operator with the kernel K n (2+x /jn 2 / 3 , Xq+ 
y/jn 2 / 3 ) defined in M.10\) is not a trace class operator in L 2 (s, oo) © L 2 (s, oo) (recall that 
A : Hi —> is a trace class operator if \\A\\\ := Tr(A* A) 1 / 2 < oo). To take care of this 
problem we follow the approach of [26] and use weighted L 2 spaces. If we take any w such that 
w^ 1 G L\ and grows at infinity not faster than exponentially, then K n is a Hilbert- Schmidt 
on L 2 (s, oo; w) ©L 2 (s, oo; w^ 1 ). The diagonal entries of K n are finite rank, hence trace class. 
Now the definition of determinant extends to Hilbert- Schmidt operator matrices T with trace 
class diagonal entries by setting 

det(J - f ) = det 2 (I - f)e' Trf , 

where Tr denotes the sum of the traces of the diagonal entries of T and the det 2 is the 
regularized 2- determinant, defined for the Hilbert- Schmidt operator T with eigenvalues t^ as 

det 2 (/-f) = l\(l-t k )e^ 

(see fl3\j . Section IV. 2). It follows from the identity for 2-determinant 

det 2 (J - fx)(I - f 2 )e TrfA = det 2 (I - f x ) det 2 (/ - f 2 ) 

that for this extended definition we still have the relation 

det(I - fx) (I - f 2 ) = det(/ - f{) det(J - f 2 ) 

(2.) Consider a rank one kernel u(x)v(y), where u € L 2 (s, oo; w 2 ) and v E L 2 (s, oo; w± 1 ) 
and u®v: L 2 (s, oo; w±) — > L 2 (s, oo; w 2 ) 

/•oo 

(u ® v h)(x) = u(x) / h(y)v(y)dy. (2-1) 



Then we have 

< 11^11^(^)11^1^2(^1) (2.2) 

Introduce the scaled kernels: 



S n (x,y) = (n 2 / 3 7 )- 1 S n (2 + x/7n 2 / 3 ,2 + y/ 7 n 2 /3) ; 

V n {x,y) = (n 2 / 3 7 )- 2 L> n (2 + x/7n 2 / 3 ,2 + y/ 7 n 2 / 3 ), 

l n {x,y) = J n (2 + x/ 7 n 2 / 3 ,2 + y/ 7 n 2 / 3 ), 

IC n (x,y) = (n 2 / 3 7 )- 1 ir n (2 + x/ 7 n 2 / 3 ,2 + 2/ / 7 n 2 / 3 ), 



(2.3) 



where S n is defined by (jl.lip and K n is defined in (jl.6p . 

We prove first assertion (ii). Observe that the determinant in (11. 14ft is the same if we 
replace the interval (2 + s/n 2//3 , 2 + e) by (s, en 2 / 3 ) and the kernel K n by 

^ x >y)-{l n ( x , y )-e(x-y) S n (y,x) J' {2A) 
Lemma 1 If we denote 

<p n (x) = n- 1 / 6 7 - 1 / 2 ^ n) (2 + xhn 2 ' 3 ), ^ n {x) = 1/2 </& (2 + x/ 7 n 2 / 3 ), (2.5) 



7 



then 

S n (x,y) = JC n (x,y) + -ip n (x)e(p n (y) +r n (x,y) 

d 1 d 

V n (x,y) = -—K, n (x,y) - -ip n (x)(p n {y) - —r n (x,y) (2.6) 
ay 2 ay 

ln(x,y) = ZK, n {x,y) + ~ei() n (x)e<p n (y) + (er n )(x,y) 

where 

\\r n (x,y)\\i, ||^r„(x,i/)||i, ||(er n )(a;,y)||i < Cn" 1 / 3 log 6 n. (2.7) 
Proof. Since (jl.36p , (jl.37p and (jl.35p imply for j < n — 2m 

n—l 

-E(^ (0 ' n) )^"V) = -</>?V) + 0(e- clog2 "), 

k=0 

we have 

2?Tl TL 1 Yl 2 TTl 

~ n E E^ n) (A)(-M ( °' n) y = E # W?V) + 0(e- c1 ^"). (2.8) 

j=0 k=0 j=0 

For n — \ <j>n — 2m we need to use the result of [21] (see Eq. (68)), according to which 
for any \p — n\ < 4 log 2 n we have 

oo 

- - (e^;M») - ^m) = n ~ l E R pd n \») + ™~%oo 

1=0 

n—l oo 

= n ^ E R vfM n) ^) + n ~ l E ^W^OO + Ti^epO*), (2-9) 

1=0 l=n 

where the remainder terms e p (fi) satisfy the bounds 

IMIl 2 [-L,L] < CrT 112 log 4 n. 
Therefore the scaled functions e p (x) = n _1//3 e p (2 + xf-yn 2 ^) admit the bounds 

11^11x2(^-1) < CrT 1 ! 2 log 4 n. 
Hence, using the definition of (fOll) and (123]) . we obtain 

n—l n—l n—l 

- n £ E^ n) (A)(^ (0 ' n) )^^ (0 ' n) ) jJt ^ n) (^)= E ^ n) wj n V) 

j=n—2m+l k=0 j'=n— 2m+l 

n-1 

+ ^ n> 0*) E ^ n) (A)(^)Ti_ 1 + ri 1 )(A^) + ri 2 )(A^), (2.10) 

j=n— 2m+l 

where (A,//) collects the terms, which appear because of the second sum in the r.h.s. of 

(2) 

(|2.9j) . and r n (A, /i) collects the remainder terms e p of (|2.9|) : 

n—l n—l oo 

j=n—2m+l p=0 l=n 
n—l n—l 

-n 2) (A,M) := E E^ (0,n) e^ n) (A)e P (M). 

j=n—2m+l p=0 
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Definition 1 We will say that some remainder kernel ?n (A, (J,) (a = 1, 2, ... ) satisfies the 
bound B with exponents k\, K2 and k^, if 

|| n -2/3 r W(2 + x / 7n 2 / 3 ) 2 + y/ 7 n 2 / 3 )||i < Cra~ Kl m K2 log K3 n. (B) 
According to the results of [7], we have 

fl -V6 7 -V2^(")(2 + a ;/ 7n 2/3) = A-(x + (n-j)/c*n 1 /3)(i + 0( ri -i/3)) ) (2.H) 
| n -i/6 7 -i/2^(»)( 2 + x / 7 n 2 / 3 )| < Ce~\ 
where c* is some constant not important for us. These asymptotic implies, in particular, that 

\\n-y^-y^ n \2 + x/ 7 n^)\\ LHw) <C, 

[|„-V6 7 - 1/2 ^ n) (2 + x/7n 2/3 )[|^( w -i) < C, (2.12) 

Using the asymptotic, (jl.32p and ()2.2|) . we obtain that rn\\, fi) satisfies (B) with k± = 1/3, 
K2 = 1, «3 = 0. 

Similarly, using (I2TTT]) . (fL32l) and ([22]), we get that ri 2) (A,/i) satisfies (B) with 

Kl = 2/3, K2 = 1, «3 = 4. 
Moreover, if we denote 

n-1 

rW(A, M ) = ne^(n) £ (^ (n) (A) - ^(A))^)-^ (2 .13) 

j=n— 2m+l 

then (USD and (|2TT0|) give us 
n— 1 n— 1 

-EE^fw (0,n) G e ^ n) (a 4 ) = K n {\,n) 

j=0 k=0 

+ ^^ n) Wt ) 1 (A)((^)- 1 e n _ 1 ,n) +r( 1 )(A,^) + r( 2 )(A, /,) + r( 3 )(A, /,), (2.14) 

where u is a vector, whose components are given by 

Ui = 1, i£ [n-2m,n-l], Uj = 0, i [n — 2m, n — 1], (2-15) 

and the remainder term r„ in view of (|2.1ip . (jl.32[) and (|2.2I) satisfies the bound (B) with 
Kl = 1/3, k 2 = 1, «3 = 0. 

Now we consider the term (see (11.36P ) 

n— 1 n— 1 

il(A >A i) = 2 £ a,4 n) (A) £ bje^in). 

k=n~2m j=n—2m 

Using (|2.9p and ()2.1ip . similarly to the above it is easy to obtain that 

A(X, ^) = \ (a, u) (b, u) ^ (A) e^ n) (jj) + r g> (A,//), (2.16) 

where u is defined in (|2.15p . and the remainder r„ satisfies the bound (B) with k\ = 1/3, 
«2 = 0, k 3 = 0. 
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Let us find (a,u)(b,u). Making transposition in (ll.36|) and taking into account that 
( - A4(°> n ))~ 1 and 2?(°' n ) are skew symmetric matrices, we get 

-(.M(M)-i = ^(^.nl^M)-!^ + }_ akb . +0 (n- 1 /2] 0gn ). 
Taking the sum of the equation with (|1.36p and applying the result to u we get 

(a, u)(b,u) = h[V^ n \(K^)- x ]u,u) = -((^ (0 ' n) )- 1 -u,P (0 ' n) u) + O(n- 1 / 2 m 2 logn). 
But it is easy to see that 

p(o,n) n = _ 6n _ 1 _|_ en _ 2m + e n - 2m -l 

Hence, 

(a,u)(b,u) = ((K^r'en-uu) - ({H^y 1 (e n . 2m + e n . 2m -i), u) + O^^m 2 log n). 

Moreover, since "P = 7£ _1 has only 2m — 2 nonzero diagonals, the standard linear algebra 
argument yields that for j < n — 2m (JZ^ ,n ^)~ l ej = Vej. Then, using (|l,30p and (|2.15p . we 
obtain 

{V{e n - 2m + e n _ 2m -i), u) = P(2) 

Finally 

(a, u)(b, u) = ((^°> n ))- 1 e n _i,M) -P(2) +0(n" 1 / 2 m 2 logn). (2.17) 

Then, combining f|2.14|> with (I2.17|) and bounds (B) for r^\\, fi) with a = 1,2,3,4, we get 
the first line of (|2.6|) . The second line of ()2.6[) can be proved similarly, if we use that (|2.1ip 
can be differentiated. To prove the last line of (12. 6|) we used that (12. lip implies that for 
\k — n\ = o{n) 

\eM^)\<Cn^ 2 . (2.18) 

Hence 

||e^(2 + x/ 7 n 2/3 )| l^^-i) < Cn- 1 ' 2 . (2.19) 

Using these bounds we get the estimates for er^' (A, fj,) and eA 2 \\,fi). The bound for 
er( 3 )(A,/i) and er( 4 )(A ;/ u) follow from ([13]) fl£2J), (|2T2|) and (fL32|) . 
□ 

Let us transform the kernel K n . We use the representation 

n— 1 oo ^oj;, 

^(a, m) = £ EE n^ (n) k* / H n) (A + ^ n V + ") + V^V + -)^ (n) (A + .)) . 

k=0j=n J ° 

(2.20) 

,9 ,9 

The representation can be obtained by taking — + — — from both sides of (I2.20P and using 

o\ on 

of (|1.35p to expand TprV'fc (A) with respect to the basis {ij)^ > }°° =1 . 

OA J 

Using the same trick as above, on the basis of (|2,lip and (|2.2p it is easy to show that 



Ti 



n—l oo 



E E v '(J {n) h* r d " (^ n) (A + ^i_ } iGu + y) 

l — n A—*, / •'O 

+^ n) (A i + ^t ) 1 (A + z,)) +ri 5 )(A,/i). (2.21) 



fc=0 j=n 
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(5) 

where (A, fx) satisfies (B) with m = 1/3, ^2 = ^3 = 0. 
Moreover, using (|1.20p . we have 

n—l 00 n— 1 00 00 

y n ■■= E E ^'(^ (n) )i^ = E E vcr )^ + o(n- x ) = E kv k + 0(0, 

fc=0 j=n k=0 j=n k=l 

where 

Vj. = — [ V'(2cosx)e ikx dx. 

On the other hand, it is evident that if we consider 

n-l 

V(x) = Vfc srn kx i 

k=0 

then 

V° = ^-V(x) +0(n- 1 ). 

ax x=0 

But it was proved in [21] (see Lemma 1) that V(x) = sinxP(2cos x). Thus we get 

V^ = P(2) + 0(n- 1 ), 

and we obtain from (|2.2ip that the kernel JC n from (|2,3p can be represented in the form 

1 f 00 

£n(a?> V) = n dz i^nix + z)(j) n (y + z) + ip n (y + z)(j) n {x + z)) + r£\x,y), (2.22) 

2 Jo 

where ||ri 6 ^||i < Cra -1 / 3 . Hence, the kernel S n is represented in the form 
1 Z" 00 

5 n (x, y) = - dz {lp n (x + z)<M2/ + Z)+ 1p n (y + 2#n(z + Z)) 

1 Jo 

+ -ip n {x)e(p n {y) + r n (x,y), \ \r n \ |i < Cmn~ 1/3 . (2.23) 

Thus we can prove that S n converges in the trace norm to Sai, repeating almost literally 
argument of [26], but using (|2.1ip instead of classical asymptotic for Hermite polynomials. 
Indeed, relations (|2.1ip yield 

lim \\<p n (- + z) - Ai{- + z)\\ L 2 {w) = lim ||V„(- + z) - Ai(- + z)\\ L 2 (w) = 0. (2.24) 

Let us prove that K n : L 2 {w) — > L 2 {w) of (|2.3p converges in the ||...||i norm to Qa* : L 2 (w) — ► 
L 2 (w). where Q A i is defined in ([L~5jl . Using (j2TL|) . (fl~9j) , and fl2T2]) we have 

ll/Cn-Q^Hx < (2.25) 

1 f°° 

2 y ^ l</,n( ' + z) " Ai( -' + z)||l2 ^ + ^ n( ' + z) _ + ^i^w) 

X {\\<Pn(- +z)\\tf(w) + [|VVi(' + z)\\l2(w) + + 2 )IIl2(«,)) dz 

<C / (||^(. + ^)-Ai(. + ^)|| WH + ||^(- + z)-AiX- + «)||ra W )e-*tfa 
J 

for some n-independent C > 0. Here we have again used ()2.1ip . implying 

\Wn{- + z)\\ L 2 {w) <Ce-\ \\i, n (. + z )\\ L2{w) <Ce- z . (2.26) 
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Now we can use the dominated convergence theorem to make the limit n — > oo in (|2.25j) . 
To pass to the limit n — > oo in ip n (x)ep n (y), remark that uniformly in y > s 

/>n 2 / 3 7 (£-2) 

eVn(y) =c Vn - / <p n (z)dz, 

1 rn 2 ^(L-2) 1/2 f L 

Vn (z)dz = — / ^(X)d\, (2.27) 



-n 2 / 3 7 (L+2) ' 2 J-L 

where L was defined in Remark [TJ But according to the results of [7J, 
n 1 / 2 / v4 n) (A)dA^0, 

J-2+n-V4 

/ f L / -2-n- 1 /4\ 

n 1 / 2 / +/ Vi n) (A)dA^0, n^oo. 

\j2+n-!/4 y 

Thus, ()2.11|) and the evenness of ip^ yield 

lim c Vn = lim n 1/2 / ^i n) (A)dA = / Ai{x)dx = 1. 

n->oo n->oo ^-n" 1 /* J_oo 

Moreover, (|2.1ip allow us to pass to the limit n — > oo in the second term of the representation 
(|2.27p of €ip n . Thus we have uniformly in y > s > — oo 

POO 

lim eip n (y) = 1—1 Ai(z)dz. 

rwoo J y 

Now implies 

lirn^ ^IC n (x,y) + ^ip n (x)e(p n (y)^ = QAi(x,y) + ^i(x) ~ J Ai(z)dz 

where the limit is understood in the || . . . ||i-norm. 

To prove that —d y K n : L 2 (w~ 1 ) — > L 2 {w) converges in the ||...||i norm to —d y QAi ■ 
L 2 (w~ l ) — > L 2 (w), we repeat the argument used in (|2.25p . taking into account (|2.2p with 
u>i = w~ l and W2 = w. Besides, we have the relations 

ll¥>n(- + 2)11^-1) <C"e- z , ||^(- + ^)|| L2(w -i) <C"e- z , (2.28) 

lim ||<p n (- + z) - + 2)||z,2 (w -i) = lim ||^ n (- + z) - Ai{- + z)\\ L 2( w -i) = 0, 

and 

lim ||^(- + z) - + 2)1^2(^-1) = lim + z) - Ai'(- + z)\\ L 2 {w -x) = 0. 

n— >oo v ' n— >oo v ' 

We obtain then 

lim (-d y IC n (x,y) - il) n (x)tp n (y)) = -d y Q A i(x,y) - -Ai(x)Ai(y), 
where the limit is understood in the ||...||i norm. 
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We are left to prove that X n : L 2 {w) — > L 2 (w 1 ) converges in the ||...||i norm to the 
operator from L 2 (w) to L 2 (u> _1 ), defined by the kernel 

J e(x - x')Q Ai (x' ,y)dx' . 

roc roc 

$ n {x) = / (p n (x')dx', ^n(x) = / ip n {x)dx. 

J X J X 



To this end denote 



Then 

1 f°° 

€(p n (x) = - / ip n (x)dx' - $ n (x) = - ® n (x), 
^ J -oo 

1 f°° 

e^n(x) = - / lj> n (x ')dx' - * n (x) = -#„(z) 
1 f°° 

eKn(x,y) = - dz(eip n (x + z)tp n (y + z) + ecp n (x + z)ip n (y + z)) 
1 Jo _ 

1 /" c 

= -- / dz(^ n (x + z)v? n (y + z) + $ n (x + z)V'„(y + 2)) + ^* n (y). 

z Jo 1 

Here the second relation follows from the fact that V^-d is an °dd function, and the third 
one follows from the first, and the second, combined with (I2.27D , Hence, repeating again the 
argument used in (|2.25p and taking into account that (|2,lip implies 

/*oo roo 

lim ||^ n (a;)— / Ai(x')dx'\\ L 2( w -i\ = lim II — / Ai(x')dx'\\ L 2/ w -i) = 0, 

n-»oo J x n->oc J x 

and 

||* n (- + 2)11^2^-1) < Ce~ z , \\® n (- + z)\\ L 2 {w -i } < Ce~ z , 

we obtain 

\xmT n {x,y) = - Q A i(x',y)dx' 

n^oo J x 

+\(^J Ai{x')dx' + J Ai(x')dx' J Ai(y')dy'^j , 

where the limit is understood in the 1 1 . . . 1 1 1 norm. Thus we have proved assertion (ii) . 

Note that we have also proved that limn^oo JC n (x, y) = QAi(x,y) uniformly in x,y £ 
(s, oo). Hence assertion (i) is also proved. □ 
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